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Spacecraft  Plasma  Modeling  Challenges 


Spacecraft  Propulsion  Relevant  Plasma: 

©  From  hall  thrusters  to  plumes  and 
fluxes  on  components 

©  Complex  reaction  physics  i.e. 
Discharge  and  Breakdown  in  FRC 

©  Relevant  Densities  often  Span 
6+  Orders  of  Magnitude 

©  Spatial  scales  of  interest  span 
/xm-100m  range 
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Spacecraft  Plasma  Modeling  Challenges 


Spacecraft  Propulsion  Relevant  Plasma: 

©  From  hall  thrusters  to  plumes  and 
fluxes  on  components 

©  Complex  reaction  physics  i.e. 
Discharge  and  Breakdown  in  FRC 

©  Relevant  Densities  often  Span 
6+  Orders  of  Magnitude 

©  Spatial  scales  of  interest  span 
/xm-100m  range 

Solution? 

Multi-Scale  and  Multi-Physics 
Adaptive  Algorithms 
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Multi-scale  Particle  Advancement 


Boris  Push: 


qAt 

v  =  v(r  -  At/2)  +  —  £(t) 
2m 

,  _  _ 

V  =  V  H - V  X  B{t) 

2m 


qAt 

2m 


1+(^|S|) 

I  qAt 

v(t+At/2)  =  v+  +  - — E(t) 

2m 

r(t-\-At )  =  r(t)  +  Atv(t  +  At/2) 


v'  X  5(0 


o  Leap-Frog  in  r  &  v 

o  Energy  Conservation  only  via 
Interpolation  in  Time 
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Boris  Push: 
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2m 

r(t-\-At )  =  r(t)  +  Atv(t  +  At/2) 


v'  X  5(0 


o  Leap-Frog  in  r  &  v 

©  Energy  Conservation  only  via 
Interpolation  in  Time 

©  Numerically  Stable  but  Drifts  and  Fails 
when  At  >  ujc/2 
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Boris  Push: 
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o  Leap-Frog  in  r  &  v 

©  Energy  Conservation  only  via 
Interpolation  in  Time 

©  Numerically  Stable  but  Drifts  and  Fails 
when  At  >  ujc/2 


Magnetic  Mirror  Test 


Multi-scale  Particle  Advancement 


Gyrokinetic  Push: 

o  Effective  for  High-B  Plasma  (i.e. 
Magnetic-Fusion) 

o  Assumes  Ac  <C  dx 

o  Looses  Phase  Information 

e  Assumes  Phase-Scatter  Diffusion 
Negligible 


Multi-scale  Particle  Advancement 


B -Frame  Push: 

o  Coordinates  Rotated  and  Aligned  to 
ZFField 
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Multi-scale  Particle  Advancement 


B -Frame  Push: 

o  Coordinates  Rotated  and  Aligned  to 
ZFField 

o  Motion  Decomposed  into  Rotation  and 
Drift 

o  Exact  Solution  in  Constant  Field 
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Multi-scale  Particle  Advancement 


B -Frame  Push: 

®  Coordinates  Rotated  and  Aligned  to 
ZFField 

®  Motion  Decomposed  into  Rotation  and 
Drift 

®  Exact  Solution  in  Constant  Field 
o  Phase  Preserved  Across  Many-cuc 
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Multi-scale  Particle  Advancement 


B -Frame  Push: 

o  Coordinates  Rotated  and  Aligned  to 
ZFField 

o  Motion  Decomposed  into  Rotation  and 
Drift 

®  Exact  Solution  in  Constant  Field 
o  Phase  Preserved  Across  Many-cuc 

What  about  Variable  Fields? 
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Multi-scale  Particle  Advancement 


B -Frame  Push  in  Variable  Fields: 

o  Consistent  Solution  Recovered  with 
Multiple  Steps 


Multi-scale  Particle  Advancement 


0 


B -Frame  Push  in  Variable  Fields: 

®  Consistent  Solution  Recovered  with 
Multiple  Steps 

o  Convergence  Order  is  Low 
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B -Frame  Push  in  Variable  Fields: 

®  Consistent  Solution  Recovered  with 
Multiple  Steps 

o  Convergence  Order  is  Low 
©  Solution  Stalls  when  dt  >  uoc 


Magnetic  Mirror  Test 
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Multi-scale  Particle  Advancement 


B -Frame  Push  in  Variable  Fields: 

o  Consistent  Solution  Recovered  with 
Multiple  Steps 

o  Convergence  Order  is  Low 

©  Solution  Stalls  when  dt  >  uoc 

©  Consider  the  Push  in  Aligned 
Coordinates... 

©  No  Mechanism  for  -A  vy  in  Step? 

(No  Bounce  in  Magnetic  Mirror?) 

©  Variable  B  and  E  via  Taylor  Series? 


Av  =  Do  •  v  +  •  E 

m 

Ax  —  Di  •  vAt  +  - - D2  •  E 

m 


Djfc  —  R  •  Ak  •  R 


Cl 

Si 
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B -Frame  Push  in  Variable  Fields: 

o  Consistent  Solution  Recovered  with 
Multiple  Steps 

o  Convergence  Order  is  Low 

©  Solution  Stalls  when  dt  >  uoc 

©  Consider  the  Push  in  Aligned 
Coordinates... 

©  No  Mechanism  for  -A  vy  in  Step? 

(No  Bounce  in  Magnetic  Mirror?) 

©  Variable  B  and  E  via  Taylor  Series? 

But  even  R  depends  on  B 
(And  is  Partially  Arbitrary!) 


Av  =  Do  •  v  +  •  E 

m 


Ax  =  Di  •  vAt  +  - - D2 

m 

Djfc  =  R  •  Ak  •  R 

> 

11 

u? 

0=^0^ 

[;]*(■?  1  3 

A0 

A, 

a/=(A,  1  »).| 
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Multi-scale  Particle  Advancement 
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B -Frame  Push  in  Variable  Fields  (cont.): 
o  Push  as  Explicit  Matrix  Operator: 

(  For  now  Ignoring  E ) 

©  What’s  Known  about  M: 

-  Reversible  through  Time  Inversion... 

-  M-1  Must  Exist 

-  M  is  a  1:1  Map  of  X0  -A-  X+ 

-  M  is  Unique  (Even  if  R  was  Not...) 

-  If  M  is  Unique,  Elements  are  Unique 

-  M  is  Constant  if  B  is  Constant 
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B -Frame  Push  in  Variable  Fields  (cont.): 
o  Push  as  Explicit  Matrix  Operator: 

(  For  now  Ignoring  E ) 

©  What’s  Known  about  M: 

-  Reversible  through  Time  Inversion... 

-  M-1  Must  Exist 

-  M  is  a  1:1  Map  of  X0  -A-  X+ 

-  M  is  Unique  (Even  if  R  was  Not...) 

-  If  M  is  Unique,  Elements  are  Unique 

-  M  is  Constant  if  B  is  Constant 

Hypothesis: 

Elements  of  M  Vary  Smoothly,  x  Nfbh  v0 

(Assuming  B  Varies  Smoothly  with  x) 
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B -Frame  Push  in  Variable  Fields  (cont.): 


o  Assuming  B  varies  Periodically  in  c j( 
and  Linearly  in  multiples  N  x  uoc...l 
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B -Frame  Push  in  Variable  Fields  (cont.): 


o  Assuming  B  varies  Periodically  in  c j( 
and  Linearly  in  multiples  N  x  uoc...l 


Multi-scale  Particle  Advancement 


B -Frame  Push  in  Variable  Fields  (cont.): 

o  Assuming  B  varies  Periodically  in  uoc 
and  Linearly  in  multiples  N  x  uoc...l 

o  Elements  of  M  take  Form: 

niij  =  my  +  +  cos(um5t)  +  sin(ujnSt)  Amj^ 
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Multi-scale  Particle  Advancement 


B -Frame  Push  in  Variable  Fields  (cont.): 


o  Assuming  B  varies  Periodically  in  ujc 
and  Linearly  in  multiples  N  x  uoc...l 

o  Elements  of  M  take  Form: 

my  =  my  +  +  cos(umSt)  Am^  +  sin(ujnSt)  Am^ 


®  Coefficients,  can  be  Sampled 
using  Original  Push  with 

=  [0, 7r/4,  37t/4,  27ra]cuc 
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B -Frame  Push  in  Variable  Fields  (cont.): 

o  Assuming  B  varies  Periodically  in  ujc 
and  Linearly  in  multiples  N  x  uoc...l 

o  Elements  of  M  take  Form: 

my  =  my  +  +  cos(um5t)  +  sin(ujnSt) 


®  Coefficients,  can  be  Sampled 
using  Original  Push  with 

At  =  [0, 7r/4,  37t/4,  27ra]cuc 

o  Comparison  of  Fit  and  Sampled  M 
from  Magnetic  Mirror 


0123456789  10 


Time,  (Gyroperiods) 


Closed  Circle  -  Original  Push  Stepping  my 
Open  Circle  -  Estimated  my  from  Sampling 
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Compounded  Push  Operator: 

o  Uniform  Field  Push  can  be  Subdivided 
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Compounded  Push  Operator: 

o  Uniform  Field  Push  can  be  Subdivided 

X(f  +  At)  =  M(A t)  •  X0 

=  m^-l\At/p)-xQ 
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Compounded  Push  Operator: 

o  Uniform  Field  Push  can  be  Subdivided 

o  Can  be  Written  as  Product  of  St 
Substeps 


X(f  +  At)  =  M(Af)  •  X0 

=  M(p-1)(A  t/p)-: 

P~  1 

=  -Xo 

k= 0 
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Compounded  Push  Operator: 

o  Uniform  Field  Push  can  be  Subdivided 

o  Can  be  Written  as  Product  of  St 
Substeps 

o  Also  True  if  varies  with  t  - 
Equivalent  to  Explicit  Stepping  with 
Smaller  At 
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Compounded  Push  Operator: 

o  Uniform  Field  Push  can  be  Subdivided 

o  Can  be  Written  as  Product  of  St 
Substeps 

o  Also  True  if  varies  with  t  - 
Equivalent  to  Explicit  Stepping  with 
Smaller  At 

o  Estimator  for  Reduces  Field 
Evaluations,  but  not  Push  Multiplier 
Operations... 
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Operation  Reduction: 

o  Given  the  Gyro-Estimated 
Compounded  Operator  M... 


~  p~l 

M  =  Yl 

k= 0 


□  >  «  g  >  <  -  ►  4  ^  ^  >0  0,0 


R.S.  Martin  (AFRL/RQRS)  Distribution  A:  Approved  for  public  release;  unlimited  distribution;  PA  #TBD.  11/27 


Multi-scale  Particle  Advancement 


0 


R.S.  Martin  (AFRL/RQRS)  Distribution  A:  Approved  for  public  release;  unlimited  distribution;  PA  #TBD.  11/27 


Multi-scale  Particle  Advancement 


0 


Operation  Reduction: 

o  Given  the  Gyro-Estimated 
Compounded  Operator  M... 

o  Can  be  Split  into  Mean  and 
Deviation 

©  The  Operator  is  Linearized 
around  M 


~  p~l 

M  =  Yl 

k=0 


p- 1 

=  Yl 

k= 0 


P~  1 

=  Yl 

k= 0 
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Operation  Reduction: 

o  Given  the  Gyro-Estimated 
Compounded  Operator  M... 

o  Can  be  Split  into  Mean  and 
Deviation 

o  The  Operator  is  Linearized 
around  M 


~  p~l 

M  =  Yl 

k= 0 
P~  1 

=  Yl  (m+amw) 

k= 0  V 
~P~  1 

=  L  Yl  (a+MAM^l)  r 

k= 0  V 
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Operation  Reduction: 

o  Given  the  Gyro-Estimated 
Compounded  Operator  M... 

o  Can  be  Split  into  Mean  and 
Deviation 

o  The  Operator  is  Linearized 
around  M 

©  And  Assuming 

> 


p- 1 


M  =  M[k 


p~  l 


=  Yl  (m+AM[A 


k= 0 

=  L 


p-2 
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Multi-scale  Particle  Advancement 


Operation  Reduction: 

o  Given  the  Gyro-Estimated 
Compounded  Operator  M... 

©  Can  be  Split  into  Mean  and 
Deviation 

a  The  Operator  is  Linearized 
around  M 

o  And  Assuming 

1 1 A||  >  ||RAMMl|| 

©  Summations  to  Series  Solution 
for  Linear  and  Oscillatory  Parts 
of  AM  — )>  Direct  Eval  of  nty 


~  p~l 

M  =  Yl 

k= 0 


p- 1 

=  Yl  (m+AM[A 


k= 0 

=  L 


p-2 
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Multi-scale  Particle  Advancement 


0 


Operation  Reduction: 

«  Given  the  Gyro-Estimated 
Compounded  Operator  M... 

o  Can  be  Split  into  Mean  and 
Deviation 

®  The  Operator  is  Linearized 
around  M 

©  And  Assuming 

1 1 A||  »  ||KAMMl|| 

©  Summations  to  Series  Solution 
for  Linear  and  Oscillatory  Parts 
of  AM  -A  Direct  Eval  of 


Magnetic  Mirror  Test 
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CR  Modeling  for  Atomic  Hydrogen 


Elementary  processes  (electron-impact  only) 

©  Collisional  excitation/deexcitation 

H(i)  +  e  //(/)+  e 

o  Collisional  ionization/recombination 

H{i)  +  e~  O  H+  +  2e~ 

9  Bound-bound  radiative  transition 

H(i)  +  hv<r-  H(j ) 
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Rate  Equations 


Rate  Equation 


P  oo 

Kj“d=Aji,  Kfl=NeJE  a  ij(£)vf(e)d, 


£ 


dNi 

dt 


E  E^'-^) 

Re Rxn  j/i 

^^K(/iVrJ,oiV, 

•/¥* 


Implicit  formulation 


dNj  AT1 


E^i-toiVi 

j+i 
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Level  Grouping 


Hydrogen 

Spectrum 
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Level  Grouping 


n,m:  bin  indices;  i,j:  level  indices 


Model 

Cons.  vars. 

Size 
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Number  density  (m  3) 


Isothermal  Heating  Test  (Te  =  3  eV) 


Comparison  of  ground  state  and  electron  number  densities 


□  S  <  1  ►  <  -1 


^  0  3  0' 


R.S.  Martin  (AFRL/RQRS)  Distribution  A:  Approved  for  public  release;  unlimited  distribution;  PA  #TBD.  16  /  27 


Number  density  (m  3) 


Isothermal  Heating  Test  (Te  =  3  eV) 


Comparison  of  excited  states  number  densities 
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Uniform  vs.  Boltzmann 
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Uniform  vs.  Boltzmann 
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Prior  Merging 


Numerous  Previous  Merge  Methods: 
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Prior  Merging 


Numerous  Previous  Merge  Methods: 


©  2:1  -  Cannot  Conserve  Energy 

(Lapenta  &  Brackbill,  JCP  1994) 


[QUES 


0 


□  >  4  gp  >  4  -  >  4  =  :  ^  '0<\CV 


LIC  RELEASE;  UNLIMITED  DISTRIBUTION;  PA  #TBD.  20  /  27 


Prior  Merging  Techniques 


Numerous  Previous  Merge  Methods: 
o  2:1  -  Cannot  Conserve  Energy 

(Lapenta  &  Brackbill,  JCP  1994) 

o  Complex  Macro-particles  with 
Internal  Energy 

(Hewett,  JCP  2003) 


Wm  =  Vm  =  'YhwiVi/Wm 

ri"0  =  (E^if"0  +  E-iV?  -  /*,„ 

Wm  (£  +  rif®))  =  E  (wiV?  +  w,-^'"'))  y 

Particle  Push  with 
Split  if  (rM)1/2  >  vm?... 

In  What  Coord-System?... 
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Prior  Merging 


Numerous  Previous  Merge  Methods: 

o  2:1  -  Cannot  Conserve  Energy 

(Lapenta  &  Brackbill,  JCP  1994) 

®  Complex  Macro-particles  with 
Internal  Energy 

(Hewett,  JCP  2003) 

o  Merge  to  Grid 

(Assous  et  al.,  JCP  2003,  Welch  et  al„  JCP  2007) 
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Conserve  Arbitrary  Moments  to  Grid 

Not  Explicitly  Conserved  Lost? 
Entropy  Generation? 

Shape  Functions? 

Grid  Dependence? 


Prior  Merging 


Numerous  Previous  Merge  Methods: 

o  2:1  -  Cannot  Conserve  Energy 

(Lapenta  &  Brackbill,  JCP  1994) 

®  Complex  Macro-particles  with 
Internal  Energy 

(Hewett,  JCP  2003) 

o  Merge  to  Grid 

(Assous  et  al.,  JCP  2003,  Welch  et  al„  JCP  2007) 

All  Introduce  Significant  Error 
and/or  Complexity 
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Conserve  Arbitrary  Moments  to  Grid 

Not  Explicitly  Conserved  Lost? 
Entropy  Generation? 

Shape  Functions? 

Grid  Dependence? 


Conservative  Merge 


Merge  to  Pair  — >>  DOF  for  Conservation: 

o  (n+2):2  yields  Exact  Mass, 
Momentum,  and  Kinetic  Energy 
Conservation 

©  Applied  Spatially  also  Shown  to 
Conserve  Electrostatic  Energy 

©  Though  Energy  Conserving, 

Still  Thermalizes  VDF 


(AFOSR  Review  2006) 
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Conservative  Merge 


Merge  to  Pair  — >>  DOF  for  Conservation: 

o  (n+2):2  yields  Exact  Mass, 
Momentum,  and  Kinetic  Energy 
Conservation 

©  Applied  Spatially  also  Shown  to 
Conserve  Electrostatic  Energy 

©  Though  Energy  Conserving, 

Still  Thermalizes  VDF 


(AFOSR  Review  2006) 

Selection  of  Near  Neighbors  in  VDF 
Limits  Thermalization 

(«  Near  Neighbor  Pairs  in  2: 1  Merges  that  Limit  Numerical  Cooling) 
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Octree  Binning 


Phase-Space  Decomposition 


o  Given  a  Set  of  Particles... 


□  >  <i  g 


^  'O  o^O 


R.S.  Martin  (AFRL/RQRS)  Distribution  A:  Approved  for  public  release;  unlimited  distribution;  PA  #TBD.  22  /  27 


Octree  Binning 


Phase-Space  Decomposition 


o  Given  a  Set  of  Particles... 
o  Particles  Binned  in  Octants 
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Octree  Binning 


Phase-Space  Decomposition 


o  Given  a  Set  of  Particles... 
o  Particles  Binned  in  Octants 
o  Octants  Recursively  Sub-Divided 
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Octree  Binning 


Phase-Space  Decomposition 

o  Given  a  Set  of  Particles... 

o  Particles  Binned  in  Octants 

o  Octants  Recursively  Sub-Divided 

o  Recursion  Halted  at  1 -Particle/Bin  or 
Other  Criteria  such  as  Bin-Density 
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Octree  Binning 


Phase-Space  Decomposition 


o  Given  a  Set  of  Particles... 

o  Particles  Binned  in  Octants 

©  Octants  Recursively  Sub-Divided 

o  Recursion  Halted  at  1 -Particle/Bin  or 
Other  Criteria  such  as  Bin-Density 

Restricts  Phase-Space  Diffusion  to 
Within  Local  Bins 
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OD-Merge  Examples 


Comparison  of  Random  vs.  Octree  Merge  Partner  Selection 
(Note:  Mass,  Momentum,  and  Kinetic  Energy  Both  Exactly  Conserved  ) 
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OD-Merge  Examples 


Comparison  of  Random  vs.  Octree  Merge  Partner  Selection 
(Note:  Mass,  Momentum,  and  Kinetic  Energy  Both  Exactly  Conserved  ) 
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Beam  in  Potential 


Collisionless  Particles  in  Well 

o  6000  Unmerged  Particles 
o  Reproduces  3-4  Orders  of  Magnitude 


Time  (ms) 


Analytical 


Unmerged 


18  20  22  18  20  22 
Number  Density,  Iog10(n  (#/m3))  Number  Density,  Iog10(n  (#/m3)) 
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Beam  in  Potential  Well 


0 


Collisionless  Particles  in  Well 

o  6000  Unmerged  Particles 
o  Reproduces  3-4  Orders  of  Magnitude 
o  Random  Merge  ->  Thermalization 
o  3000  First  Point,  1500  First  Cross 
o  Bi-Maxwellian  Specifically  Difficult 
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Beam  in  Potential  Well 


Collisionless  Particles  in  Well 

2 

((I)) 

2 

A 

o  6000  Unmerged  Particles 

1.8 

(l) 

1.8 

A 

o  Reproduces  3-4  Orders  of  Magnitude 

1.4 

1.4 

x 

o  Random  Merge  ->  Thermalization 

I1'2 

A 

I1'2 

V 

o  3000  First  Point,  1500  First  Cross 

1  i 

0.8 

A 

1  1 

0.8 

A 

o  Bi-Maxwellian  Specifically  Difficult 

0.6 

x 

0.6 

w 

o  Octree  Merge  Significantly  Better 

0.4 

0.2 
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-0.5  0  0.5 
Position  (m) 
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Number  Density,  logl  0(n  (#/m3)) 
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Beam  in  Potential  Well 


Collisionless  Particles  in  Well 

©  6000  Unmerged  Particles 
o  Reproduces  3-4  Orders  of  Magnitude 
o  Random  Merge  ->  Thermalization 
©  3000  First  Point,  1500  First  Cross 
©  Bi-Maxwellian  Specifically  Difficult 
©  Octree  Merge  Significantly  Better 
©  Merge  &  Split  Adapts  Particle  Count 
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Beam  in  Potential  Well 


Collisionless  Particles  in  Well 

o  6000  Unmerged  Particles 
q  Reproduces  3-4  Orders  of  Magnitude 
o  Random  Merge  ->  Thermalization 
o  3000  First  Point,  1500  First  Cross 
o  Bi-Maxwellian  Specifically  Difficult 
o  Octree  Merge  Significantly  Better 
o  Merge  &  Split  Adapts  Particle  Count 

o  Computational  Particles  per  Cell 
Vastly  Different 


Unmerged 
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DC-Breakdown 


DC-Diode  Test  Case: 

o  Full  3D  Electrostatic-PIC 
®  Averaged  to  ID  XT-Plot 
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Position,  (mm)  Position,  (mm) 


Number  Density,  (#/m!) 


Number  Density,  (#/m!) 


Control  Merged 


□  >  <i  g 


FOR  PUBLIC  RELEASE;  UNLIMITED  DISTRIBUTION;  PA  #TBD.  25  /  27 


DC-Breakdown 


DC-Diode  Test  Case: 

o  Full  3D  Electrostatic-PIC 
o  Averaged  to  ID  XT-Plot 
©  250V  Cathode  — >>  Anode 
©  MCC-Ionization  Collisions 
©  Secondary  Emission  at  Cathode 
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Time,  (ns) 


Number  Density,  (#/m3)  Number  Density,  (#/m3) 

Control  Merged 
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DC-Breakdown 


DC-Diode  Test  Case: 

o  Full  3D  Electrostatic-PIC 
o  Averaged  to  ID  XT-Plot 
©  250V  Cathode  — >>  Anode 
©  MCC-Ionization  Collisions 
©  Secondary  Emission  at  Cathode 
©  Chain-Branching  Needs  Merge 
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DC-Breakdown 


DC-Diode  Test  Case: 

o  Full  3D  Electrostatic-PIC 

o  Averaged  to  ID  XT-Plot 

©  250V  Cathode  — >>  Anode 

©  MCC-Ionization  Collisions 

©  Secondary  Emission  at  Cathode 

©  Chain-Branching  Needs  Merge 

©  Merge  Overhead  Rapidly 
Negligible 
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DC-Breakdown 


DC-Diode  Test  Case: 

o  Full  3D  Electrostatic-PIC 

o  Averaged  to  ID  XT-Plot 

©  250V  Cathode  — >>  Anode 

©  MCC-Ionization  Collisions 

©  Secondary  Emission  at  Cathode 

©  Chain-Branching  Needs  Merge 

©  Merge  Overhead  Rapidly 
Negligible 

©  Control:  Parts/Cell  oc  Density 
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Time,  (ns) 


Computational  Particles  per  Cell 

Control  Merged 
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DC-Breakdown 


DC-Diode  Test  Case: 

o  Full  3D  Electrostatic-PIC 

o  Averaged  to  ID  XT-Plot 

©  250V  Cathode  — >>  Anode 

©  MCC-Ionization  Collisions 

©  Secondary  Emission  at  Cathode 

©  Chain-Branching  Needs  Merge 

©  Merge  Overhead  Rapidly 
Negligible 

©  Control:  Parts/Cell  oc  Density 

o  Merge:  Parts/Cell  Much 
Reduced 
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DC-Breakdown 


DC-Diode  Test  Case: 

o  Full  3D  Electrostatic-PIC 

o  Averaged  to  ID  XT-Plot 

©  250V  Cathode  — >>  Anode 

©  MCC-Ionization  Collisions 

©  Secondary  Emission  at  Cathode 

©  Chain-Branching  Needs  Merge 

©  Merge  Overhead  Rapidly 
Negligible 

©  Control:  Parts/Cell  oc  Density 

©  Merge:  Parts/Cell  Much 
Reduced 

©  Despite  Identical  Densities 
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Control  Merged 
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AFRL/RQRS  M&S  Future  Work 


0 


Integrate  R&D  w/  Production  TODO: 

®  High-Order  Fluid/MHD  GPU  Models 

(Le/Cole*/Bilyeu  PhD  Research) 

o  GPU  Accelerated  Chemical  Kinetics  / 
CR  Ar-Ne-Xe-Molecular  Models 

(Le/Cole*/Kapper*  PhD  Research) 

o  Phase- Space  Reconstruction/Vlasov 

(Martin/B  ily  eu/TB  D) 

o  Implicit  /  Multiscale  GPU- Accelerated 
PIC 

(Lederman/Gimelsheins/Martin/TBD) 


*Note:  Former  Co-op  Student  Work  to  be  Integrated  into  Framework 
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End 

Thank  You 

Questions? 
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